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Abstract 

We calculate the black hole entropy in Loop Quantum Gravity as a 
function of the horizon area and provide the exact formula for the 
leading and sub-leading terms. By comparison with the Bekenstein- 
Hawking formula we uniquely fix the value of the 'quantum of area' 
in the theory. 

^ . The Bekenstein-Hawking formula gives the leading term in the entropy 

lO ! of the black hole in the form 

°- <? 1A m 



o 



where lp is the Planck length and A is the black hole horizon area. According 
to the common belief, the entropy is always connected with the logarithm of 



the number of microscopic states realizing a given macroscopic state. The 
vL| ! fact that the entropy in (1) is proportional to the area (and not as is usual to 

the volume) has led to the formulation of the so called holographic principle. 
It was however difficult to find the microscopic states that could account for 
^ ■ such an entropy. 

The problem was attacked in two different approaches - in string theory 
[l]-[5] and in Loop Quantum Gravity [6]- [10] (see also the review [11]). The 
latter approach is based on a quantum theory of geometry. The basic ge- 
ometric operators were introduced in [12, 13] and a detailed description of 
the quantum horizon geometry was introduced using the U(l) Chern-Simons 
theory in [10] where one can find the detailed discussion of the states on the 
horizon of the black hole (see also [14]). However, the procedure introduced 
in [10] for state counting contained a spurious constraint on admissible se- 
quences and the number of the relevant horizon states is underestimated. 
The correct method of counting was proposed in [14] but that analysis pro- 
vides only lower and upper bounds on the number of states. The purpose 
of this paper is to rectify this situation. We will provide the correct value 
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of the Barbero-Immirzi parameter that is needed to obtain agreement with 
the Hawking-Bekenstein formula for large black holes. Furthermore, since we 
are able to calculate the number of states to sufficient accuracy, we also rig- 
orously obtain the precise sub-leading quantum correction to the Hawking- 
Bekenstein formula. The physical aspects of this result will be discussed 
elsewhere [15]. 

As it was shown in [14] the proper counting of states is given by all 
sequences (of arbitrary length) of m 8 6 Z/2, m 8 ^ such that 

\/KI(KI + 1) < « (2) 

i 

where 

° — 7, — w (3) 

and 7 is a parameter introduced in [16]. Addditionally we have to impose 
the condition 

^m i = 0. (4) 

i 

It is our task to calculate the number of sequences satisfying (2) and (4). 

We start with the sequences without (4) imposed. If we denote by N(a) 
the number of sequences satisfying (2) then we can split the counting into 
the first number m x and the rest and write the recurrence relation 

N(a) = 6{a - y/l/2) (2N(a - y/3/2) + 2N(a - V2) + . . . 

+2N(a - VKKKI + 1) + • • • + 2 \\fAa 2 + 1 - ll) (5) 



where the symbol [. . .] denotes the integer part. The first term on the RHS 
corresponds to sequences with at least two elements and m 1 = ±1/2, the 
second to sequences with at least two elements and m\ — ±1 and so on and 
the last term to sequences consisting of just m 1 and nothing else. 

We now assume (we will later prove that it is indeed the case) that the 
leading part of the solution of the recurrence relation (5) is for large a given 
by 

N(a) = Ce 2niMa (6) 
where 7m is a constant to be determined. 
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Plugging the form (6) into (5) we get 



y^2e~ 27r7M v /fc(fc+2)/4 (7) 



k=i 



where we neglected terms vanishing in the limit a — > oo. The solution to this 
equation can be found numerically and reads 

7m = 0.23753295796592 ... (8) 

We see that the coefficient of growth 7m is between the lower and upper 
bounds given in [14] i.e. ^ < 7m < It seems that 7m defined by 

(7) cannot be calculated analytically or given as a solution of any "simpler" 
equation. In the present paper we will also use 7 where 

7m := 2tt 7 = 1.492463591462379 ... (9) 

Since the condition (4) (as we will later prove) gives only power-like cor- 
rections to N(a) therefore for large a we have from (6) 

S = lnN(a) = ^^ + 0(\nA). (10) 

The comparison with the Bekenstein-Hawking formula (1) gives therefore 
the unambiguous result 

7 = 7m (11) 

so this physical requirement uniquely fixes the value of 7 and hence the 
'quantum of area' in the theory. 

One may note that the exponential ansatz (6) works in many other cases 
for example if the number of states for a given spin j was (2j + 1) instead of 
2 one would have in (7) (k + 1) instead of 2 in front of the exponential and 
the appropriate 7m would be equal to 0.273985635 . . .. 

We now prove that the assumption (6) is actually correct and at the end 
we add the condition (4) - the final result will not modify (11) and will enable 
us to find the sub-leading (logarithmic) correction to the formula for entropy 
(10). 

We introduce the Laplace transform of N(a): 

PCX) 

P(s) : = / daN(a)e- sa . (12) 
Jo 
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The transform is well defined since we know that N(a) is piecewise continuous 
and of exponential order i.e. it is bounded by Ke~^ a for some K and (3 ([14]). 
Integrating the relation (5) over a from to oo we get an exact result 



9 V 00 - Sv /fc(fc+2)/4 

P{s) = . ZZ - k = l6 . . (13) 

^l_ 2 ^ =1 e-v / ») 

If all the poles of the Laplace transform are distinct and of finite order then 
we can make analytical continuation to the complex half-plane Re(s) > 0, 
excluding positions of the poles, i.e. beyond the region allowed by the def- 
inition (12). We now use the well known fact that the growth of N(a) is 
determined by the poles of P(s) i.e. 

N(a) = res ^e Sia + 0{a n ) (14) 

Si,Re(sj)>0 

(we assumed here that all the poles are simple - otherwise we would have to 
add some powers of a in front of the exponentials). Therefore our problem 
boils down to determination of such complex Sj that satisfy 

G(si) -\ = (15) 

where we defined 

oo 

G{ S ) ; = ^ e -V fc (fc+D/4. ( 16 ) 
k=l 

It is important to analyze the distribution of zeroes of (G(s) — |) in more 
detail. Since G(s) is complex analytic in the open domain Re(s) > all the 
zeroes are distinct, of finite order and there is no point of accumulation of 
zeroes inside this domain. Therefore we have to analyze the domain's border 
i.e. the imaginary line Re(s) = 0. To analyze the limit when Re(s) \ 
define a new function G r (s) by subtracting two analytic (for Re(s) > 0) 
functions 

G r (s) := G(s) - e- s(fc+1)/2 " \ E ( 17 ) 

k=l k=l 

The function G r (s) is well defined in the whole domain Re(s) > 0, continuous 
and doesn't have any singularities in this domain. On the other hand two 
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subtracted functions have a well defined limit when Re(s) \ so we can 
write 



G(s) = G r (s) + 



o-s/2 



se 



-8/2 



■\n{l-e- s ' 2 ). 



(18) 



e s/2 _ I 4 

We see that apart from singularities at s = 4irin, the function G(s) has 
a well defined limit when Re(s) \ and in the domain Re(s) > it is 
complex analytic so in the whole domain Re(s) > there can be no points 
of accumulation of zeroes - therefore all zeroes of the function G(s) — \ in 
the whole domain Re(s) > must be distinct and of finite order. 
The only real zero is given by (9): 



s = 7m- 

Close to this zero the function P(s) behaves as 

Cm 



P(s) 



s -7m 



where 



C 



M 



0.509202564 



2 7 mG"(7m) 

Therefore the leading behaviour for large a is given by 



(19) 



(20) 



(21) 



N(a) = C M e 



7MO 



(22) 



as was to be shown. 

The complex zeroes of the function G(s) — \ for Re(s) > 1. and |Im(s) 
100 are listed below: 



< 



Re(s) 


Im(s) 


1.49246359. . . 





1.22393017... 


± 22.1530069... 


1.41016352... 


± 35.9362749. . . 


1.30363023. . . 


± 58.0472739. . . 


1.18587535... 


± 71.8281531... 


1.17654302... 


± 79.9673974. . . 


1.07191106... 


± 87.4120423. . . 


1.21660746. . . 


± 93.9713059... 
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It is easy to prove that all the complex zeroes (i.e. poles of P(s)) with nonva- 
nishing imaginary part must have smaller real part than 'Jm so in comparison 
to the leading behaviour (22) they give exponentially small (and rapidly os- 
cillating) contribution to N(a). 

To illustrate the method in the cases where the number of sequences is 
explicitly known let us give two examples. For sequences of half-integers 
with the condition \rrii\ < a the number of sequences is 3t 2a l - 1 and by 
the method described above we get 

P M = 7(^h) < 23 > 

and indeed the real pole is equal to 2 In 3. For sequences of half-integers with 
the condition XXI m i| + 1/2) < a the number of sequences is |(2l 2a ] +1 ± 1) — 1 
(where + is for even and — for odd [2a]) and 

*M = (24) 

and indeed the positive real pole is equal to 2 In 2. The actual numbers 
of sequences in both cases indeed show that the leading behaviour is given 
precisely by the largest positive real pole of P(s). 

It can be useful to discuss how different ways of defining the sequences 
are reflected in the properties of P(s). We have essentially three different 
"statistics" . 

• Boltzmann statistics (used in this paper): all sequences satisfying the 
condition are allowed and are treated as different. Therefore P(s) is 
given by eq. (13): 

P(s) = . Z1 - k = l6 7 (25) 

*(l-2Er=ie-^ fc(fc+2)/4 ) 

which has the largest positive real pole equal to 7m- 

• the intermediate statistics (between Boltzmann and Bose): it requires 
that > \rrii\ but treats for example (+|, — |, +|, . . .) as different 
from (— |, +|, +|, . . .) - then we have 

oo 1 

P(s) = TT 7 N (26) 

t\ (l - 2e- s v / ^)7i) 
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which has the largest positive real pole given by k — 1: s = -^=-. 

• Bose statistics: it requires not only that > \rrii\ but treats for 

example (+|, — |, +|, . . .) as the same as (— |, +|, +|, . . .) i.e. the 
allowed sequences have the form (|, . . . , — |, . . . , +1, . . . , —1, . . .). Then 
we have 

= II 7 ; ^ ( 27 ) 

t x (l - e -V*(*+ 2 )/ 4 ) 

and this P(s) does not have a pole with Re(s) > at all so iV(a) must 
have milder growth than e ca . However when s — > the function P(s) 
behaves as P(s) ~ where ((s) is the Riemann's zeta function. 

Hence we have accumulation of zeroes at s = and the inverse Laplace 
transform gives the exponential growth with \fa: 



7V(a)~a c eV / ^ (28) 

(this type of statistics and therefore similar formula gives for example 
the degeneracy of states in string theory). 

At the end we impose the condition (4) i.e. J^m^ = 0. We introduce 
N(a,p) as the number of sequences satisfying 

a/KKK + 1) < a (29) 

i 

and _ 

S ^ j m i =p (30) 

i 

where p G Z/2 and a > a/N(M + !)■ We can write the recurrence relation 

as 



N(a,p) = e{a-V3/2)6(a-y/\p\(\p\ + l)) x 

(V(a - p - 1/2) + (iV(a - v^/2, p + 1/2)+ 

+N(a - V2,p - 1) + N(a - V2,p+1) + . . . + 1) . (31) 



where the first term inside the brackets on the RHS corresponds to sequences 
with at least two elements and mi = +|, the second to sequences with at 



least two elements and mi = — | and so on and the last term to sequences 
consisting of just mi (which must be equal to p). 

Next we define P(s,uj) as the Fourier transform with respect to p and 
Laplace transform with respect to a of N(a,p): 

POO 

P(s,u):=Y] / daN(a,p)e iujp e~ sa . (32) 
„^^ /n Jo 



Then summing and integrating eq. (31) we get 
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P ( S) u) = - s J2 eW*(*+2)/4 l-2£ e W*(*+2)/4 C os(^/2) (33) 

k=i V fe=i / 

where we neglected some contributions from the endpoints {p ~ a) what will 
be justified a posteriori. Therefore we have to find the zeroes of the function 
inside the bracket on the RHS of (33). With oo = we know that the only 
real zero (i.e. the one that gives the leading behaviour) is equal to s = 7^. 
Expanding around s = 7m we get 

oo 

1 - 2 e-W*(*+2)/4 (l - ( s - j M )y/k(k + 2)/A- A; V/8 + 0(w 4 )) = 0. 
fc=i 

(34) 

Hence the pole of P(s, cu) is given by 

s = 1m - (3 moo 2 + 0(u 4 ) (35) 

where 

Pm = -i^L___ = 0.475841255 . . . (36) 

Therefore the leading behaviour of N(a,p) is given by 

N(a, P) = ^r J duJ e ~ iuJP e^ Ma -? MUj2a = e 2 ^ Ma e - p2/ ^ Ma \ (37) 

As we see the endpoints p ~ a are exponentially small justifying (33). Also we 
see a posteriori that the terms 0(oj a ) from (35) would give in (37) corrections 
of relative size 1/a so for large a we can neglect them. 
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Note that the gaussian distribution with respect to p in (37) with the 
mean square deviation of the order of y/a can be expected by the analogy 
with the random walk where the average distance is zero and the mean square 
deviation is of the order of square root of the number of steps. 

Imposing the condition p = we get the final result 

N(a, 0) = ° M e 2 ™. (38) 



Therefore the entropy is given by 



7m A HA/P P ) 
47 I] 



5 = ln^(a,0) = ^^-^™ + O(l) (39) 



so for large A/lp we confirm (10) as the leading behaviour and we can un- 
ambiguously both get the physical value of 7 equal to 7^ (given in (8)) 
and resolve the controversy (summarized in [17]) about the coefficient of the 
logarithmic term fixing it to — |. 
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